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CHEMICAL SPECIES SUPPRESSION FOR MRI IMAGING USING 
SPIRAL TRAJECTORIES WITH OFF-RESONANCE CORRECTION 

BACKGROUND OF THE INVENTION 
[0001] The present invention relates to a method of fat suppression in MRI imaging 

and more particularly to a method of fat suppression using Dixon techniques with spiral 
trajectories. 

[0002] OfiF-resonance effects (e.g., field inhomogeneity, susceptibility, chemical shift) 

cause artifacts in magnetic resonance imaging (MRI). The artifacts appear as positional shifts 
along the readout direction in rectilinearly sampled acquisitions. Usually, they are 
insignificant, because of short readout, times, in normal spin echo (SE) and gradient echo 
(GRE) sequences. However, off-resonance artifacts sometimes appear as severe geometric 
distortion because of the relatively long readout time in echo planar imaging (EPI). 
[0003] Over the past decade, spiral imaging techniques have gained in popularity due 

to their short scan time and insensitivity to flow artifacts. However, off-resonance effects 
(e.g., ~3.5ppm shift between fat and water protons) cause blurring artifacts in the 
reconstructed image. Most spiral off-resonance correction methods proposed to date are 
difficult to apply to correct for blurring artifacts due to the fat signals since the fat- water shift 
is typically much greater than that due to BO inhomogeneities across the field of view (FOV). 
As such, off resonance artifacts remain one of the main disadvantages of spiral imaging. 
[0004] Currently, off-resonance artifacts are most commonly avoided by use of 

spatially and spectrally selective radio-frequency (RF) excitation pulses (SPSP pulses) since 
they excite only water spins, thereby eliminating the off-resonance fat signal and hence avoid 
artifact generation. Yet, SPSP pulses usually take a significant fraction of the TR of a spiral 
sequence, particularly when a sharp transition band is designed between two chemical shift 
frequencies. The time duration of SPSP pulses is sometimes even equivEdent to the spiral 
readout duration. This can be considered as one of the current limitations to reducing scan 
time in spiral imaging, and thus, alternative approaches to fat-water separation in spiral 
imaging are required. 

[0005] Excitation of only water spins could be achieved through application of 

chemical shift presaturation pulses (e.g., CHESS pulses) prior to normal spatially selective 



excitation. However, the effectiveness of these frequency selective RF excitation pulses is 
dependent on the main magnetic field (Bo) inhomogeneity. 

[0006] Altematively, Dixon techniques have primarily been investigated for water-fat 

decomposition in rectilinear sampling schemes. In the original Dixon technique, water and fat 
images were generated by either addition or subtraction of the "in-phase" and "out-of-phase" 
data sets. Water and fat separation is unequivocal using this technique when magnetic field 
inhomogeneity is negligible over the scanned object. However, when Bo inhomogeneity 
cannot be neglected, the original Dixon technique fails to accurately decompose water and fat 
signals. Therefore, modified Dixon techniques using three data sets (i.e., three-point Dixon 
(3PD) technique) or four data sets were developed to correct for Bo inhomogeneity ofiF- 
resonance effects and microscopic susceptibility dephasing. New versions of the Dixon 
technique use two data sets with Bo inhomogeneity off-resonance correction (i.e., two-point 
Dixon (2PD) technique). The water-fat decomposition performance is almost equivalent to 
that of the 3PD technique although off-resonance frequency estimation of this technique 
would be imstable for voxels with nearly equal water and fat signal intensities. The 
advantage of these multiple-point Dixon techniques over spectrally excited RF pulses is that 
water-fat separation can be achieved even in the presence of Bo inhomogeneity. This 
advantage is of notable importance because neither tissue-induced local magnetic field 
mhomogeneity nor extemal applied magnetic field inhomogeneity can be completely 
removed by shimming. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0007] The invention may take form in certain components and structures, preferred 

embodiments of which will be illustrated in the accompanying drawings wherein: 
[0008] Fig. 1. is a simplified sequence diagram of the spiral three-point Dixon (spiral 

3PD) technique; 

[0009] Fig. 2. illustrates an example of Diocar J^? plots that have two local minima; 

[0010] Fig. 3. illustrates a selection process of the correct local off-resonance 

fi-equency using a region growing £dgorithm; 

[0011] Fig. 4. is a simplified sequence diagram of the spiral two-point Dixon (spiral 

2PD) technique; 

[00 1 2] Fig. 5. illustrates water and fat magnetization vectors in the transverse plane at 

the onset of the spiral gradient in the spiral 2PD technique; and 
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[0013] Fig, 6. illustrates variable density spiral (VDS) trajectories for Dixon 

techniques with one acquisition. 

DETAILED DESCRIPTION OF THE INVENTION 
[00 1 7] It is to be understood that the specific devices and processes illustrated in the 

attached drawings, and described in the following specification are simply exemplary 
embodiments of the inventive concepts defined in the appended claims. Hence, specific 
dimensions and other physical characteristics relating to the embodiments disclosed herein 
are not to be considered as limiting. 

[00 1 8] The long readout time of spiral trajectories leads to off-resonance signals that 

blur into neighboring pixels; spins firom multiple ofiF-resonance fi-equencies can all contribute 
to the voxel signal. However, it can generally be assumed that Bo inhomogeneity is smoothly 
varying across the FOV (this assumption will be referred to as 'assumption (i)'). Thus, the 
average off-resonance firequency in any pixel is usually close to the true local Bo field 
strength. This concept is exploited in the conventional method to create an off-resonance 
frequency map in spiral imaging, in which the phase difference is taken between two images 
(with different TE's) even though both images are already blurred by off-resonance effects. 
[0019] In addition to assumption (i), it is also assumed that: (ii) Only two chemical 

shift species, i.e. water and fat, are considered, and their spectra are both sufficiently narrow 
so that little spectral overlap occxirs, and, (iii) Signal intensity differences due to T2* decay 
among data sets with different TE's are negligible. While water and fat are used as the two 
chemical species in the examples provided herein, it should be appreciated that the invention 
can be applicable to any other two suitable chemical species. 

Spiral 3PD technique 

[0020] Suppose that three k-space data sets with TE differences (ATE (s)) are 

acquired using a normal spatially selective RF pulse for each TR (Fig.l). As shown, the TE 
difference between the first and the second acquisitions and that between the second and the 
third acquisitions are set to equal (ATE). 

[002 1 ] The signals at each pixel in the reconstructed images (iSo, Si and Sz) firom these 

data sets can be expressed as: 

5'o=PF'+F' [1] 
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5, = (JV'+F' exp(/^^^ )) expm [2] 
= iW'+F' exp(i2^^„ )) exp(/2(^) [3] 
where W represent water signals blurred by local Bo inhomogeneity off-resonance frequency 
/(Hz), F' are fat signals blurred by local Bo inhomogeneity and chemical-shift ofif-resonance 
frequencies /+ f/ai (Hz), and ^ and ^j-^ are the phase shifts due to Bo inhomogeneity and 
chemical-shift off-resonance effects during ATE. That is, 

^ = 2n:f ATE, W 
respectively. 

[0022] A quadratic equation with respect to expOV) can be derived from Eqs.[l-3] to 

provide an estimate of the frequency field map: 

5, (expOV,^ ) + 1) ± yjs,^ (exp(/V^„ ) + !)' - ^S^S^ expQV^ ) 
expOV) = — — — ; LOJ 

In practice, it is difficult to select one appropriate root from the two roots given by Eq. [6]. 
Moreover, at the pixels where image SNR is low, both roots of Eq.[6] may be inaccurate. 
[0023] Alternatively, several different predetermined frequencies _f} (Hz) for a 

sufficiently wide frequency range are substituted into Eqs.[l-3], and W and F' are solved 
for each frequency j^. The calculated W and F' (Eqs.[7,8]) are compared with the observed 
signals, So, Si and 52. The objective is to find the jfj that provides the best fit of calculated W 
and F' to the observed signals. Given a frequencyj^, W and F' can be solved from any two 
equations of Eqs.[l-3]. To make signal contributions from So, Si and ^2 equal to the 
calculated W and F' , they are computed as follows: 

f'j=^{f,,j+f;,j+f^j) m 

_ Sq+S.+S:,- F'j (1 + exp(/(^^. + <t>fa, )) + exp(/2((^^. + (^^^ ))) 
^ ~ 1 + expOVy ) -H exp(/2(*y ) 

where 

^. ^exp(i2<^,)5o-exp(zV,)5. (F' derived from Eqs.[l, 2]), 

exp(/2(^,)(l-exp(/if>^„,)) 

p' = exp(/V,).y,-53 (F- derived from Eqs.[2, 3]), 

''■^ exp(/2j*^ ) exp(/9}^ )(1 - exp(iV/« )) 
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P' expii2^j)S,-S, (F* derived from Eqs.[3.1]). 

'"•^ exp(i2(!J^)(l-exp(i2$>^„)) 

At each pixel, variation between the observed and calculated signals, Dpi^i (Eq.[9]) will be 
minimized aty = / if the frequency// is the closest to the actual Bo inhomogeneity off- 
resonance frequency out of all the test frequencies^, 
^pixw = I ^0 - (W'j+F'j ) l+l 5, - (W'j+F'j exp(zV^„))exp(iV^) I 
+ 1 S2 - (W'j +F'j expii2^f„ )) exp(z 2<*^ ) | [9] 
In practice, Dpuei may not take the minimum value wheny = / at those pixels where image 
SNR is low. To reduce the errors from noise, a small multi-pixel window is set about each 
location and the summation ofDpucei within the window is measured: 
D - YD tlO] 

^local " pixel 
window 

[0024] Even if noise effects are completely eliminated, a plot of Z>wwith respect to 

fj sometimes takes two local minima as shown in Fig.2. This can be understood since Eq.[6] 
gives two roots. Since both roots in Eq.[6] satisfy Eqs.[l-3]. Di„cai may not be the actual local 
field strength. To select a correct frequency A a region growing algorithm is applied. The 
region growing algorithm was originally proposed as a phase unwrapping algorithm. In our 
method, the initial frequency map for the region growing algorithm can be determined as the 
value of fj that minimizes Diocai where Diocai takes the single minimum in the Diocai-fj plot. 
This defmes a 'frequency determined' region. 

[0025] Alternatively, a specific ATE can be chosen for simplification. The TE of the 

fu^ data set wUl be an integer multiple of ATE. For example, a ATE = 2.2 ms for 1 .5 T can 

be chosen and: 

5, =(jr'-F')exp(/<^) 

^2 =(fr*+F')exp(/2(!)) 

52/5o =exp(/2(*) 

[0026] The second step is to find the correct frequency// at each pixel in a 'frequency 

to-be-determined region' which abuts the 'frequency determined' region found in previous 
steps. Figure 3 illusfrates the process used to select a correct local frequency/. The shadowed 
region denotes the 'frequency determmed' region. At each pixel that borders on the shadowed 
region, the value of/ chosen is that which creates a local minima in the Aocar/plot, and, is 
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the closest to the average local frequency of the neighboring pixels in the frequency 
determined region. In Fig.3, suppose that the frequency of the pixel with a circle inside is to 
be determined and that Diocai of this pixel takes more than one local minima The value of^ at 
which Diocai takes local minima and is the closest to the average local value in the nearby 
frequency determined region is selected asfi for the pixel. The frequency determined region is 
thereby expanded so that the frequency field map can be created for the entire scanned object 
region. 

[0027] After the frequency field map is obtained, the deblurred water and fat images 

FT and F are reconstructed based on the frequencies indicated in the frequency field map at 
each pixel location. Specifically, the water and fat images decomposed with the frequency^?, 
i.e., W\ and F\ derived from Eqs. [7, 8], are demodulated with demodulation frequencies^/ 
and // + ffat to create the locally deblurred images W, and F,, respectively. The entirely 
deblurred water image W can be reconstructed by combining the deblurred regions of water 
images W, from each local frequency, A in the frequency field map. The entirely deblurred 
fat image F can be reconstructed in a similar fashion. 

[0028] To simplify the method described above, ATE can be chosen as r, the time 

during which fat spins precess by 180° out of phase with respect to water spins. This 
condition provides for a more efBcient determination of the frequency map, With this 
condition, the signals at each pixel in the reconstructed images {Sq, Si and S2) from these data 
sets can be expressed as: 

5o = W'+F' t^^l 
5, =(lfF'-F')exp(/(z5) pa] 
= (PT'+F') exp(i2^) [3a] 
where W represent water signals blurred by local Bo inhomogeneity off-resonance frequency 
/(Hz), F' are fat signals blurred by local Bo inhomogeneity and chemical-shift off-resonance 
frequencies /+ f/a, (Hz), and ^ is the phase shift due to Bo inhomogeneity off-resonance 
effects during ATE. (The term f/a, can also be referred to as/,*,/, or/ for chemical species 
other than fat) That is: 
^ = 1 n f bTE = l7if x. 

[0029] Local off-resonance frequencies are determined. 14 is obtained as: 

li, = Arg{SJS^). [5a] 
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[0030] Even though signals blur into neighboring voxels due to the local off- 

resonance effects in the spiral image, from assumption (i), 24 in Eq.[5a] gives the phase shift 
due to the true local Bo inhomogeneity off-resonance frequency / at each pixel. Phase 
unwrapping is then performed to obtain the correct (/> at each pixel using the region growing 
method with a manually selected seed point as described above, /can be successively 
determined from Eq. [4a] . 

[0031] Water-fat signal decomposition is performed at each pixel based on the 

determined frequency map. From Eqs.[la-3a], one obtains: 

W = ((So + exp(-/2<^)) 12 + S, exp(-/V)) / 2 , [6] 
F'= ((So + -^2 exp(-/2<i})) / 2 - S, exp(-i^J)) / 2 . [7] 
[0032] De-blurring needs to be performed on the W and F' images. This is achieved 

through the use of the frequency-segmented off-resonance correction method. This method 
reconstructs several images using different demodulation frequencies and selects the most de- 
blurred region from the stack of reconstructed images under the guidance of a frequency map. 
The demodulation frequencies used to de-blur W are those given by the frequency map. The 
demodulation frequencies used to de-blur F' are the sum of the chemical-shift off-resonance 
frequency and the local frequencies given by the frequency map. The use of multiple 
demodulation frequencies enables de-blurred water and fat images, and F to be 
reconstructed. 

Spiral 2PD technique 

[0033] Although the spiral 3PD technique can achieve water-fat separation with off- 

resonance deblurring, the scan time is prolonged since three data sets are acquired. Here, a 
newly developed 2PD technique for spiral trajectories is described. 
[0034] hi the spiral 2PD technique, as shown in Fig.4, the TE's of the first and second 

data acquisitions are nrand (n+l)r, respectively, where « is a positive integer and ris the 
time during which fat spins precess by 180° out of phase with respect to water spins. 
[0035] The signals at each pixel in the reconstructed images (So and S 0 from these 

data sets can be expressed as: 

So =Pf^'+F', t^^^ 
S, =(PJ^'-F')expOV), t^2] 
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where the definitions of JF' and F' are the same as in the spiral 3PD technique, and ^ is the 
phase shift due to Bo inhomogeneity ofif-resonance frequency /(Hz) during r(s). That is, 

Eq. [12] can be obtained by simply substituting expOV^^ ) = -1 in Eq. [2]. 
[0036] Since it can be assumed that an RF pulse tip angle into the transverse plane is 

the same for both water and fat spins in each voxel, water and fat magnetization vectors in 
the transverse plane should be aligned with each other at the onset of the spiral gradient when 
TE is an even-integer multiple of rand they are in the opposed phase when TE is an odd- 
integer multiple of r (Fig.5). In other words, if FF' and F' are deblurred by k-space data 
demodulation with the correct local off-resonance frequencies (the water and fat images 
obtained this way are defined as W and F , respectively), the orientation of the two vectors 
W and F will be identical when TE is an even-integer multiple of r, and they will be 
opposite when TE is an odd-integer multiple of r. 

[0037] In Fig. 5, the long and short arrows represent the water and fat magnetization 

vectors, respectively. is the RF penetration angle. ^ is the precession angle for the water 
magnetization during TE. (a) When TE is an even-integer multiple of r, the water and fat 
magnetization vectors are aligned with each other at the onset of the spiral, (b) When TE is 
an odd-integer multiple of r, the water and fat magnetization vectors are in the opposed phase 
at the onset of the spiral. 

[0038] For water-fat signal decomposition, as was done in the spiral 3PD technique, 

several predetermined demodulation frequencies jfj (Hz) are substituted into Eqs. [1 1,12]. W 
and F' are solved for each demodulation frequency fj (defined as W/ and F/, 

respectively.): 

^ Soexp{i(^j) + S, [14] 
^ ~ 2exp(i>^) 

^. ^ Soexp(i<pj)-Si [15] 
^~ 2exp(i>j.) 

where = 27zf, • t . 

[0039] The deblurred water and fat images Wj and F,. can be obtained after k-space 

data demodulation with the demodulation frequencies^ and^ + f/ai for W/ and F/, 
respectively. To determine the correct Bo inhomogeneity off-resonance frequency at each 
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pixel, the orientations of two vectors Wj and Fj are compared. When n is even (odd), the 
value of ^ that makes two vectors Wj and Fj aligned (opposed) to each other is selected as 
the correct local off-resonance frequency (the correct frequency is defined as /J)- As is 
evident, this process of local off-resonance frequency determination simultaneously 
reconstructs the final deblurred water and fat images, W and F . 

[0040] In practice, since some voxels contain predominantly either only water or fat 

spins, the vector alignment property described above is difficult to use. Hence, the following 
quantity is measured: 



The absolute value of this quantity is minimized when the two vectors Wj and Fj are 
aligned/opposed or when either Wj or F,. predominantly exists. The similar concept was 
taken advantage of in an earlier rectilinear 2PD technique. 

[0041] As explained m the spiral 3PD technique, it is difficult to determine the 

correct frequency in image regions with low SNR. This is also true with the spiral 2PD 
technique. Therefore, in our algorithm, we first select the pixels of which the absolute values 
of do not take sufficiently small values at any 7. It is considered that the SNR of these 
pixels are low. Their frequencies are not evaluated in the following algorithm but are 
interpolated afterward based on the determined frequencies in their vicinity. 
[0042] After the 'frequencies to be evaluated' regions are determined, a similar 

algorithm to the spiral 3PD technique is used to reduce the errors in frequency evaluation: 
The summation of exp(j<Dp within a small window centered on each pixel is measured and 
the real part of this quantity is extracted: 

^/»c./=Re(Ze'^P('^>)) ^^^^ 

window 

Note that the summation in Eq. [17] is a complex sum. is expected to be maximized 
wheny = /. However, since often shows periodic-like patterns withy, even in the 
sufficiently high SNR regions, a plot of -fj often has. more than one peak with their 
magnitudes close to one another. In our algorithm, all the values of ^ at which Pw takes 
local maxima in the Piocar fjV^oX, are first chosen at every pixel (these^ are defined asfp). 
Then the region growing algorithm is performed to select the correct frequency/^ in a manner 
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similar to the spiral 3PD technique. In other words, of all the values of fp, the one selected is 
that which is the closest to the average local frequency of the neighboring pixels in the 
frequency determined region. Unlike the spiral 3PD technique, the region growing algorithm 
starts with a manually selected single seed point in the spiral 2PD technique. 

Algorithms for 2PD and 3PD techniques with multi-coil data sets 
[0043] Multi-element surface coils are often used to obtain images with higher SNR 

than would be achieved with a single larger coil. Each surface coil usually has a small region 
of signal sensitivity. Thus the image reconstructed from each set of individual coil data, 
shows non-uniform signal intensity over an FOV. The Bo inhomogeneity frequency map 
derived from an individual coil may not be accurate for the region where the image SNR is 
low. This subsection describes how to apply the 3PD and 2PD spiral reconstruction 
algorithms described previously when data are acquired from multi-element surface coils. 
Diocai in Eq. [ 1 0] is redefined as 

m 

to extend the spiral 3PD algorithm from a single data set to one for multi-coil acquisitions. 
Here, n is the total number of coils and Diocafim is D,ocai defined in Eq. [10] of the m-th coil 
data. The frequency field map can be determined by applying the same algorithm described in 
the previous subsection of the spiral 3PD technique using Aoca/as defined in Eq. [18]. 
[0044] For the spiral 2PD technique, at each demodulation frequency j^, the quantity 

expressed as Eq. [16] is calculated for each coil data. They are combined at each pixel, 
weighted by the signal intensities of the reconstructed images. Hence, <D^. is redefined as: 



= ArgiY, i\WjC„ I + \FjC„ I) • exp(/ • 2Arg 



)) [19] 



where n is the total number of coils, WjC„ and F,.c„ are the water and fat images 
reconstructed from the w-th coil data with the predetermined demodulation frequency^. Pw 
is computed from Oj defmed in Eq. [19] using the same algorithm as explained in the 
previous section of the spiral 2PD technique. The method to determme the frequency field 
map from the newly defined Piocai is also the same. 

One-acquisition spiral Dixon techniques using variable density spiral trajectories 
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[0045] Ideally, only one acquisition would be needed for fat/watw decomposition. 

An efficient spiral off-resonance correction metiiod with only one acquisition has been 
proposed. This metiiod is called 'off-resonance correction using variable density spirals 
(ORC-VDS)'. In tiiis method, odd- and even-numbered spiral interleaves have slightiy 
different TE's and tiie central portion of k-space is oversampled using variable density spiral 
trajectories. A Bo inhomogeneity field map can be calculated by taking tiie phase difference 
between tiie two low resolution images reconstructed from tiie data of odd- and even- 
numbered spiral interleaves. 

[0046] The extensions of tins metiiod to tiie spiral 3PD and 2PD techniques are quite 

straightforward. Figure 6 shows tiie spiral trajectories tiiat oversample tiie inner regions of k- 
space witii tiiree times (a) and twice (b) higher sampling densities tiian the outer parts. 
[0047] In Fig.6a, tiie solid, dashed, and dotted spiral trajectories have tiiree different 

TE's witii tiieir TE differences constant. In Fig.6b, the solid and dashed spiral trajectories 
have different TE's witii tiieir TE difference equal to r. hi botii cases, off-resonance 
frequency field maps can be derived from tiie low spatial frequency data using tiie same 
algoritiims as tiie spiral 3PD and 2PD techniques described above. It can be expected tiiat 
tiiese spiral acquisition schemes make it possible to perform tiie algoritiims of spiral 3PD and 
2PD techniques witii only one acquisition. The Dixon techniques taking advantage of tiie 
sampling schemes shown in Fig.6a and b will be referred to as 'variable density spiral 3PD 
(VDS-3PD) technique' and 'variable density spiral 2PD (VDS-2PD) technique' from tiiis 
point on. A description of tiie reconstruction algoritiims of VDS-3PD/2PD techniques 
follows. 

[0048] hi botii tiie VDS-3PD and VDS-2PD techniques, a frequency field map is 

derived from only low spatial frequency data using the spiral 3PD and 2PD algorithms, 
respectively. Water-fat decomposition and k-space data demodulation are performed for the 
low spatial frequency data based on tiie frequencies indicated in tiie frequency field map. 
High spatial frequency data are also demodulated witii tiie demodulation frequency indicated 
in tiie frequency field map and will be.added to tiie demodulated low spatial frequency data. 
[0049] In the VDS-3PD technique, tiie high frequency data of three different TE's 

(tiiey are defmed as Sho, Shi and Shi) are combined so tiiat tiiefr phases are consistent witii 
one anotiier. hi otiier words, when tiie demodulation frequency is fi (Hz), tiie high spatial 
frequency data of tiie tiiree TE's to be added to tiie water image, are combined as (it is 
defmed as Shw): 
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Shw = 5^0 + ^''i txp(-i27f,^TE) + Sh, exp(-/2#, • , [20] 

and those that will be added to the fat image are combined as (it is defined as Shf): 
Shf = SK + Sfh exp(-/2;r(/, + //« )^TE) + Sh, exp(-/2;r(/, + ^, ) • 2ATE) . [21] 
[0050] Similarly, in the yDS-2PD technique, the high frequency data of two TE's 

{Sho and Sh{) to be added to the water image are combined as: 
Shw = Sho + Sh^ exp(-i27f,T) , t^^^ 
and those to be added to the fat image are combined as: 

Shf = Sho - Sh^ exp(-i27f,T). t^^l 
[0051] In both VDS-3PD/2PD techniques, Shw and Shf are demodulated with the 

demodulation frequencies /, and/,+^„ respectively. These demodulated high frequency data 
are added to the low frequency water and fat images that are akeady demodulated by the 
same frequencies yj and^+ ffa„ respectively. As described above,^„, can be referred to as/^j^ 
or^ for applicable chemical species other than fat. 

[0052] Both the spiral 3PD and 2PD techniques were implemented for in-vivo 

imaging experiments. All experiments were performed using a 1.5-Tesla Siemens Sonata 
scanner (Siemens Medical Solutions, Erlangen, Germany). In these experiments, axial brain 
and pelvis images were acquired from a healthy volunteer using a quadrature head coil and 
four-element phased array surface coils, respectively. All procedures were done under an 
institutional review board approved protocol for volunteer scanning. 
[0053] The following sequence parameters were the s^e for all the spiral sequences 

used in the experiments: for the brain image experiments, there were 20 spiral interleaves. 
FOV 240x240mm, slice thickness 10mm, flip angle 13°, spiral readout time 16ms. and TR 
25ms. For the pelvis imaging experiments, there were 20 spiral interleaves, FOV 
390x390mm slice thickness 1 0mm, flip angle 1 3°, spiral readout time 1 5ms, and TR 25ms. 
[0054] For the spiral 3PD technique, TE' s were set to 3 .0(1 st)/4.5(2nd)/6.0ms(3rd) in 

both thebrain and pelvis imaging experiments. ATE wassetto 1.5ms as the maximum SNR 
is attained when = 120°. as described in ref.l8. 

[0055] For the spiral 2PD technique, TE's were setto 2.2(1 st)/4.4ms(2nd) in both the 

brain and pelvis imaging experiments. 

[0056] The normal spiral sequences with SPSP pulses were also implemented for 

comparison. 1-4-6-4-1 binomial pulses, were used for excitation. The total flip angle for on- 
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resonance spins was 16". Two acquisitions were performed for off-resonance correction. TR 
was 33nis and TE's were set to 6.0/7.5(ms) in both the brain and pelvis imaging experiments. 
[0057] The VDS-3PD technique was implemented for the axial pelvis image. The 

image was scanned from the same volunteer using four-element surface coils. The sequence 
parameters were: 18 interleaved spirals, TE 3. 0/4.5/6. 0ms (each TE was shared by 6 
interleaves), TR 25ms, flip angle 1 3^ FOV 390mm x 390mm, slice thickness 10mm, and the 
radius of the over-sampled region was 40% of kmax- 

[0058] For reconstruction of the brain images, k-space data re-gridding was 

performed to a Cartesian grid. The modified Block Uniform Resampling (BURS) algorithm 
was used for k-space re-gridding. For reconstruction of the pelvis images, the matrix 
lescaling algorithm (the equivalent algorithm was also proposed as 'next-neighbor re- 
gridding') was used to facilitate reconstruction of the multiple k-space data sets. For both 
water and fat images, the image reconstructed from each coil data was combined using the 
sum-of-squared method. 

[0059] The predetermined demodulation frequencies are ranged from -200Hz to 

+200HZ with the frequency resolution 1 OHz (i.e., 4 1 demodulation frequencies in total) in all 
image reconstructions. The window sizes to compute were set to 9x9 pixels and 5x5 
pixels for the brain and pelvis image reconstructions, respectively. The window sizes to 
compute Piocai were set to 9x9 pixels and 5x5 pixels for the brain and pelvis image 
reconstructions, respectively. 

[0060] Both spiral 3PD and 2PD techniques successfully perform off-resonance 

correction with water-fat signal decomposition. As mentioned, one of the main advantages of 
the Dixon technique is that uniform fat suppression can be achieved across an FOV in the 
presence of Bo inhomogeneity. Both the spiral 3PD and 2PD techniques lead not only to 
uniform fat suppression but they avoid unwanted water signal suppression when Bo 
inhomogeneity exists. 

[0061] Undesirable water signal suppression could be reduced by the. use of SPSP 

pulses with a sharper transition band between the water and fat frequencies. However, such 
SPSP pulses usually take a relatively long time for excitation. As mentioned, the time 
duration of SPSP pulses is sometimes equivalent to or may be longer than die spiral readout 
time. Therefore, the minimum possible TE and TR are determined by the lengths of the SPSP 
pulses. On the other hand, since normal small flip angle RF pulses are used and hence short 
TE's can be set in both the spiral 3PD and 2PD techniques, these techniques can provide a 
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way to reduce the scan time in spiral imaging. 

[0062] A five-pxilse binomial scheme was used for spatial-spectral excitation in 

experiments. It was about 9ms in duration. The length of a SPSP pulse in other previous 
studies using spiral trajectories was about 1 6ms. The spiral 3PD and 2PD techniques employ 
normal small flip angle RF pulses and the excitation usually takes only 2-3ms for each TR. If 
short TE's (2~4ms) are used in the newly proposed techniques, TR can be possibly reduced 
by approximately 10ms from the spiral acquisitions with SPSP pulses. For example, total 
acquisition time may be reduced by about 200ms (or almost half) in 20 interleaved spiral 
trajectories. The newly proposed techniques can be used to improve temporal resolution in 
some applications of spiral imaging, such as BOLD unaging and dynamic breast imaging. 
[0063] As the spiral 2PD technique achieves almost the same performance as the 

spiral 3PD technique, it is suggested that the spiral 2PD technique be used if a shorter total 
acquisition time is required. However, in addition to better image SNR, the spiral 3PD 
technique has a couple of advantages over the spiral 2PD technique. When a frequency field 
map is determined, no Fast Fourier Transform (FFT) needs to be performed in the spiral 3PD 
technique; this can be understood from Eq. [9]. On the other hand, a number of FFTs are 
required in the spiral 2PD technique. Specifically, at each predetermined demodulation 
frequency fj, two FFTs (to demodulate W/ and F,') are needed and hence the total number 
of FFT required is twice the total number of predetermined frequencies. As mentioned, some 
of these FFTs may be uimecessary to reconstruct the final deblurred . images W and F . 
Hence, the spiral 3PD technique is computationally more eflFicient than the spiral 2PD 
technique since the FFT is the most computationally intensive part in our algorithm. Another 
advantage of the spiral 3PD technique over the spiral 2PD technique is that the initial 
frequency field map can be determined without manual seed point selection when the 
frequency field map is determined using the region growing algorithm. As explained above, 
the initial frequency field map consists of the pixels of which Aoca/-^ plots have single 
minimum. Alternately, the spiral 2PD technique needs manual seed point selection to start the 
region growing algorithm when a frequency field map is determined as is also the case with 
the rectilinear 2PD technique. 

YDS Dixon techniques 
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[0064] In the conventional spiral acquisition method with SPSP pulses, two data sets 

with different TE's usually need to be acquired to correct for ofif-resonance blurring artifacts. 
The necessity of two data sets diminishes the advantage of spiral imaging as a fast acquisition 
method. When non-negligible motion is involved between the first scan and the subsequent 
scan, the off-resonance correction algorithm may fail in spiral imaging because of motion- 
dependent misregistration between the two images. Similarly, accurate water-fat 
decomposition may not be achieved in the Dixon technique if there is motion-dependent 
misregistration among the reconstructed images with different TE's. VDS-3PD/2PD 
techniques have overcome these drawbacks. 

[0065] As is evident in the algorithms of VDS-3PD/2PD techniques described above, 

since high fi-equency data can not be separated into water and fat signals, when the combined 
high fi-equency data of water Shw (Eqs. [20, 22]) are added to the low resolution water image, 
fat signals in the high fi-equency data are also added to the low resolution water image. 
However, as evident from the difference between Eqs. [20 and 21] and the difference 
between Eqs. [22 and 23], the phases of fat signals in the high frequency data of different 
TE's are not consistent with one another when water signals in the high frequency data are 
phase-consistent. Moreover, when the water image is reconstructed, the fat signals in the high 
frequency data are demodulated not by their demodulation frequency yj + but by the 
demodulation frequency of water signalsyj. Therefore, the high firequency components of the 
fat signals are smeared out in the water image. As the total signal amounts of the high 
fi^quency components are usually quite small compared with those of the low frequency 
components, the artifacts due to high frequency components of the fat signals are 
inconsequential in the water image. Similarly, since water signals in the high frequency data 
are spread out in the fat image, there are usually no apparent artifacts in the fat image caused 
by the high spatial frequency water signals. 

[0066] The spiral 3PD and 2PD techniques proposed here can achieve not only water- 

fat signal decomposition but also ofif-resonance blurring artifact correction. The SPSP pulses 
commonly used in spiral imaging may not lead to uniform fat suppression and they may also 
yield undesirable water signal suppression in the presence of Bo inhomogeneity. Furthermore, 
they usually require a long time for excitation. These disadvantages of SPSP pulses can be 
overcome by the techniques we have proposed which provides excellent performance as well 
as reduced acquisition time. Although they require additional computational costs to the 
conventional spiral off-resonance correction method, faster data acquisition is possible than 
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in conventional spiral acquisitions that use SPSP pulses. One-acquisition spiral Dixon 
techniques using variable density spiral trajectories have also been shown to be simple 
extensions of the spiral 3PD and 2PD techniques, thereby avoiding the need for multiple 
separate acquisitions. The spiral Dixon techniques proposed here offer new and alternative 
approaches to achieve effective fat signal suppression with reduced scan time. 
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